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ABSTRACT 

The early B supergiant LMC star BI 108 is photometrically variable with a unique light curve; 
two strong periods are present in an almost precise 3:2 resonance. We collected spectroscopic 
data at VLT/UVES, sampling the supercycle of 10.733 days in ten epochs. We find spectral 
signatures for a SB2 system consisting of two massive BI supergiants orbiting at the orbital 
period of 5.366 days. The shorter periodicity resembles the light curve of an eclipsing binary 
with periodicity 3.578 days that is not detected in the data. We discuss possible causes for the 
short periodicity and conclude that the quadruple system is the more plausible hypothesis. 
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1 INTRODUCTION 

The object LMC SC9-127519 (OGLE 05 1343. 14-69 1837) was dis- 
covered as a variable star with unique photometric properties dur- 
ing our multiperiodicity search in the OGLE-II database ( ISzyman^ 
|ski|2005^ . Independently this object was announced as two binary 
systems in 3:2 resonance by |Ofir|f2008| >. With the coordinates, the 
magnitude (V = 13.3) and the color (B - V = -0.142), the star 
was identified as BI 108 in the list of Brunei et al. (1975) , for wh ich 
they give a spectral type of bl-j^ |Schmidt-Kaler et al. ( 1999 1 list 
bl:II: for the same star. 

BI 108, in fact a binary or multiple system rather than a sin- 
gle star, is likely one of the brightest members of a young open 
cluster; NGC 1881 (logt < 6.7, |Pietrzynski & UdalskipOOO) . Al- 
though the resonance is not absolutely perfect, such a near resonant 
and long-term stable period relation is very rare. In order to inves- 
tigate the photometric behavior of BI 108 all available photometric 
data were used, as well as a set of 20 UVES observations to deter- 
mine the spectroscopic behavior. Brief reports of our advances in 
the understanding of BI 108 were presented in recent conferences 
( |Kolaczkowski et al.|20T0l|Rivinius et al.|201ll . 



2 SPECTROSCOPIC OBSERVATIONS AND ARCHIVAL 
DATA 

We secured 20 spectra of BI 108 between October 1, 2008 and lan- 
uary 3, 2009, with the ESO VLT/UVES. The instrument was used 
in the 437/760 setting, covering most important He! and Balmer 
lines and several other transitions of particular interest for hot 



* E-mail: rmennick@astro-udec.cl. Based on ESO Proposal 382.D-0311 
^ Spectral classes are written in lower-case letters when they have been 
determined from multicolor photometry. 



stars. At a total of ten epochs two subexposures were taken with 
tcxp ~ 1480 sec each, with a slit width of 0.8 arcsec, giving an ef- 
fective resolving power A/ A A of about 60000. The resulting typi- 
cal S/N is between 90 and 100 for both arms of the instrument. 

The UVES standard pipeline ran smoothly on our spectra pro- 
viding wavelength calibrated 1-D spectra virtually free of instru- 
mental artifacts. The spectra were continuum normalized and he- 
liocentric corrected. The observing log is given in Table[T] 

Apart from OGLE-II, data were extracted as well from 
the MACHO (MACHO object ID 79.5378.25) and the OGLE- 
III database ( |Udalski et aT] [20081 LMC 111.2.62038). The 
MACHO data span seven and a half years (HID 2 448 826. 14 
to 2451541.96), the OGLE-II data almost four years 
(HID 2 450 457.76 to 2451872.75), and the OGLE-III data 
over almost eight years (HID 2 452 167.64 to 2454951.52). 

The photometric data of all three databases were joined for the 
analysis. Since the bands are slightly different, this was done the 
following way: for the time series analysis presented in Sect. |3.I| 
the mean magnitude was subtracted from each dataset. For the 
MACHO data, only observations where simultaneous blue and 
red points exist were accepted. The two measurements were aver- 
aged to reduce the individual uncertainty. Note that this combined 
dataset was only used to determine the periods, not for any other 
analysis. For the actual disentangling of the light curve we chose 
the OGLE-II data set only, in order to avoid problems due to the 
different filter bands of each dataset. 



3 PHOTOMETRIC PROPERTIES OF THE SYSTEM 

3.1 Time series analysis and disentangling 

The time series analysis returns two very strong periods, given in 
Table|2] For reasons that will be shown later, the longer one will be 
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Figure 1. Upper panel: The OGLE-II I-band light curves for both periods of the BI 108 system (upper panel). For each period the contribution of the respective 
other one has been removed. Lower panel: The original light curve folded with the supeiperiod. In all cases HJDq = 2 451 163.8915. Dashed lines indicate 
phases of the spectral observations. 



called the orbital period (P2) and the shorter one the resonant pe- 
riod (Pi). From the light curve alone, however, this is by no means 
obvious. For none of the two periods is the associated variation si- 
nusoidal in shape, for both the minima are more pronounced and 
sharper than the maximum associated variations. In fact both vari- 
ations are of double wave nature. 

The stronger variation is associated with the shorter period, 
Pi = 3.577958 d, with a peak-to-peak amplitude of about 0.18 
mag. From the shape alone it is tempting to attribute it to an eclips- 
ing binary nature (Fig.[T| lower curve in the upper panel), however, 
we will show later that at least spectroscopy shows no sign of any 
binary with that period. 

The longer period shows weaker variations, with about 0.05 
mag peak-to-peak, and has the value of P2 = 5.366520 d (Fig.[T] 
upper curve in the upper panel). In the case of P2 the two half waves 
are indistinguishable in photometry, however, and the double wave 
nature is shown by spectroscopy. 

The most interesting feature of both periods is that they have 
extremely similar superperiods, i.e. they are almost fully reso- 
nant with a ratio of 3 : 2, or 2 x P2 = 3 x Pi = 10.733d. 
Using the law of propagation of uncorrelated errors we find 
P2/P1 = 1.499883 ± 0.000006. This differs from 3/2 = 1.500000 
by 0.000117 ± 0.000006, i.e. at the 117/6 = 19.5 sigma level. The 
difference between the two superperiods is significant but so tiny, 
that we will come back to it in the discussion, but ignore it for 
the rest of this work. Instead, we use our OGLE-II superperiod of 
10.73309 d ( (Rivinius et al.pofl ) to analyze variability with the su- 
percycle. The combined light curve phased with that superperiod 
is shown in the lower panel of Fig. [T] The minima of the com- 
bined curve are governed by the stronger variation with Pi, while 
P2 only modulates these minima, i.e. making them more or less 



Table 1. Summary of spectroscopic observations discussed in this paper. 
We adopted HJD=2451 163.8915 as epoch for the deepest minimum. All 
spectra are 1480 s exposures, the dates refer to the beginning of the exposure 



Spectrum 


HJD 


superphase 


phase (Pi ) 


phase (P2) 


(label) 


(-2450000) 


10.73 d 


3.58 d 


5.37 d 


BIl 08-01 


4 834.63285 


0.0031 


0.9318 


0.0077 


BI108-02 


4 834.65057 


0.0047 


0.9368 


0.0110 


BI108-03 


4742.76051 


0.4434 


0.2546 


0.8882 


BIl 08-04 


4742.77823 


0.4450 


0.2595 


0.8915 


BIl 08-05 


4746.83997 


0.8234 


0.3947 


0.6483 


BI108-06 


4746.85768 


0.8251 


0.3997 


0.6516 


BI108-07 


4 803.60272 


0.1120 


0.2593 


0.2255 


BI108-08 


4 803.62044 


0.1137 


0.2642 


0.2288 


BI108-09 


4761.66650 


0.2048 


0.5386 


0.4111 


BI108-10 


4761.68422 


0.2065 


0.5435 


0.4144 


BI108-11 


4796.80764 


0.4789 


0.3601 


0.9594 


BI108-12 


4796.82536 


0.4806 


0.3651 


0.9627 


BI108-13 


4744.75145 


0.6289 


0.8110 


0.2592 


BI108-14 


4744.76917 


0.6305 


0.8160 


0.2625 


BI108-15 


4741.79717 


0.3536 


0.9853 


0.7087 


BI108-16 


4741.81489 


0.3553 


0.9902 


0.7120 


BI108-17 


4745.76182 


0.7230 


0.0934 


0.4475 


BI108-18 


4745.77953 


0.7246 


0.0983 


0.4508 


BI108-19 


4747.72260 


0.9057 


0.6414 


0.8128 


BI108-20 


4747.74031 


0.9073 


0.6464 


0.8161 



deep. We choose the date of a deepest minimum as epoch, so that 
the ephemeris for the superperiod is 

Pminiigtit .super = HJD 2 451 163.8915 + Ex 10.73309 

and the ephemerides for orbital and resonant period are 
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Table 2. Results of the time series analysis on the combined photometric 
dataset over 16 years. Amplitudes are MACHO b, MACHO r, and OGLE-II 
I-band, respectively. The two superperiods, 3 X Pi and 2 X P2 , differ by a 
tiny, yet significant fraction. 



period 
[d] 



superperiod 
[d] 



amplitude 
[mag] 



Pi 3.577958(14) 
P2 5.366520(8) 



10.73388(4) 
10.733040(16) 



0.172, 0.171,0.178 
0.050, 0.047, 0.054 



Table 3. Results to the fit of the combined spectrum at $2= 0.81. 



Parameter 


star#l 


star #2 


Tefi [K] 


27 500 ± 250 


27 500 ± 250 


log 9 


3.75 ± 0.05 


3.5 ± 0.05 


V sini [kms~^] 


116±5 


116±5 


■"rad [kmS-1] 


401 ± 2 


75 ± 2 


relat. flux contrib. 


1.0 ±0.02 


0.88 ± 0.02 



Tminiight,2 = HJD2 451 163.8915 + Ex 5.366520 
and 

Tminiight.i = HJD2 451 163.8915 + Ex 3.577958 
respectively. 

3.2 Stellar parameters 



Table 4. System properties of BI 108 derived from the spectroscopic data. 



Using E{B - V) ^ 0.156 (for NGC 1881, |Pietrzynsk i & Udalski] 
|2000[> and the interstellar extinction law derived by[Schlegel et al. 
\l99% we calcu lated {B - V)o = -0.30 (indicating bOII/III, 
|Fitzgerald|1970) , {V - K)o = -0.95, (J - K)o = -0.28 and 
{H — K)o = —0.18 mag. The corresponding figures for a bl lb 
star are —0.19 mag, —0.55, —0.09 and —0.03 mag, respectively 
( |Fitzgerald| 1970[ IKoomneef] 1983^ . Although variabihty could af- 
fect some of these measurements (though not the MACHO colors, 
which are taken simultaneously), it is difficult to avoid the conclu- 
sion that the infrared colors are bluer than those of a typical bl (or 
even bO or late o type) supergiant. 

If we use V = 13.3 for the whole system, and assume two 
equal stellar components, we get V = 14.05 for a single star. Us- 
ing the distance modulus for the LMC of 18.5 mag jPietrzyhsE] 
[etaD'2009), and the NGC 1881 reddening E{B - V) = 0.156 
(Pietrzyhski & Udalski 2000), we obtain Mv = -4.96. Using 
a bolometric correction of BC = —2.70 for a star with T^s — 
28 000 K, log g = 3.0 ,t;turb = lOkms"^ and Z/Zq = 0.5 
( |Lanz & Hubeny||2007[ >, as representative for BI 108, we get an 
absolute bolometric magnitude of Mboi = —7.66. From the evolu- 
tionary tracks for Z = 0.004 stars by |Maeder & Meynet|p001| l this 
luminosity implies a total mass for the system of about 40 Mq . The 
same exercise with four stars gives a total mass of about 60 Mq. 

The MACHO color term, Bm — Rm, does not show any peri- 
odic signature, except a signal with a period of 23h56m, i.e. clearly 
related to Earth's rotation, and apart from this no other variability 
at all above the statistical error, which is 0.01 mag. 

Photometric data were retrieved from 2MASS and Spitzer ob- 
servational databases and converted to fluxes using the virtual ob- 
servatory tool VOS^^( Bayo et al.| 2008), and compared to pub- 
lished spectral energy distributions. There is no evidence for any IR 
excess in BI 108, i.e. a significant flux contribution from a cooler 
component in the visual range is firmly excluded. 



4 SPECTROSCOPIC VARIATIONS 

4.1 Spectral classification of the components 

In the spectra the system obviously has SB2 nature and contains 
two similar components with spectral features of B 1 Ib-II stars. Al- 



Property 


Value 


Property 


value 


1 
e 

a sin i 


0.75 
0.08 
93 
41.5 P0 


Ml sin^ i 
M2 sin^ i 
Msin^ i 


19 Mq 
14.5 Mq 

33.5 Mq 



though both spectra are virtually identical in terms of lines present 
and V sin i, in the combined light one component produces, on av- 
erage, slightly weaker lines (and as will be seen from the RV curve 
is as well less massive). In the following, this star as being the sec- 
ondary of the current configuration, will be denoted star #2, while 
the other one will be star#l. 

At quadratures the spectral lines are almost separated and can 
easily be compared to synthetic and other template spectra. A com- 
parison with the Bl Ib-II galactic star HD54764 gives an almost 
perfect match between line profiles (even in v sin i = 125kms~^). 

The line ratio Hel 4387/011 4416 is a useful luminosity indi- 
cator[^The average value is 0.99 ± 0.09 for star#l and 1.1 ± 0.3 
for star- #2, in agreement with a luminosity class I/II derived from 
the visual classification above. 

In spite of what is said above about one star on average being 
brighter than the other, in some of the quadrature phases the two 
spectra appear of almost equal strength, i.e. the relative contribu- 
tions from each star to the equivalent width (EW) are variable. This 
behavior is discussed in detail below, here it is sufficient to say that 
the spectral types, understood as defined by spectral line ratios, for 
instance Hel 447 1/MgII 4481, do not change significantly with the 
variable EW. 

In order to obtain individual stellar parameters for both com- 
ponents, the combined spectru m was modeled using the grid of 
synthetic spectra published by Lanz & Hubeny (2003 2007 1^ 
A multi-parameter fit of those models was done to the UVES 
data taken at good separation of both components (HJD 2 454 747, 
02 = 0.81). The fit parameters were effective temperature, sur- 
face gravity (in the range TcS 25 kK - 30 kK, log g 3 - 4), radial 
velocities and v sin i for each component, as well as the relative 
flux contribution of secondary component, fixed as unity for the 
primary (see Table[3|. The calculations were done for models with 
the average LMC metallicities and microturbulence velocities 2 and 
lOkms^^. While the hydrogen and helium lines are rather insen- 
sitive, the higher microturbulence gives the better fit for the metal 
lines (Fig.|2]l. 



^ http://www.laeff.inta.es/svo/theory/vosa/ 



^ http://nedwww.ipac.caltech.edu/level5/Gray/frames.html 

^ http://nova.astro.umd.edu/Tlusty2002/tlusty-frames-models.html 
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Figure 2. The combined spectram at <I>2= 0.81 along with the best fit whose 
parameters are given in Table 3. 




Phase (orbital period) 

Figure 3. H/3 hne (o) and Sim 4553, 4568, 4574 lines (A) radial velocities 
obtained after Fourier spectral decomposition (filled symbols star#l, open 
symbols star#2). The two outliers in the secondary's Sim RVs occur when 
the spectra are of somewhat lower quality in that region. In other cases, 
when only one type of symbol is visible this is due to agreement resulting 
in the two symbols to be plotted above each other. The entire curve is shifted 
by the 7 velocity (227 it 3 kms~^) and centered around km s^^. 



4.2 Time series analysis and disentangling 

We measured RVs for both components by finding line centers 
with deblending algorithms, as well by cross correlation of selected 
spectral regions with templates, but by far the best result was ob- 
tained by Fourier disentangling the spectra wi th the virtual obser - 
vatory program KORELp1( jHadravajl995(|Skoda & Hadrava|2010 1. 
Although robust in spectral disentangling and binary parameter de- 
termination, KOREL does not provide the zero point of velocity 
neither meaningful errors for the system parameters. 

The corresponding velocities show that P2 — 5.366 d is 
clearly a binary orbital period and suggests a near circular orbit 
(Fig. |3j. The period Prcs = 3.577 d is actually not present in the 
RVs. The fit for H/3 RVs gives Ki = 169 ± 10 kms"^ and 
K2 — 224 ± 10 kms^^, i.e. a mass ratio q of about 0.75. 

The analysis gives a small non-zero eccentricity of about 
e — 0.08 with a longitude for the periastron of cj = 93 degree. 



Similar figures were obtained by using the strongest Sim lines, 
providing extra confidence on the above results. While this small 
eccentricity is compatible with e — 0.00 plus non-zero noise, we 
note the system is relatively young, so a not fully completed circu- 
larization is entirely possible. Residual radial velocities, corrected 
for orbital motion with these parameters, are not sorted with Pi, 
i.e. that period does not reflect in the orbital dynamics of the binary 
system. The entire set of derived parameters is given in Table|4] 

We calculated independent system parameters with PHOEBE 
l |Prsa et al.|20TT| using the disentangled orbital light curve and as- 
suming Ti = T2 = 27500 K, P = 5.36654 d, g = 0.75, e = and 
no third light. Also we assumed a detached configuration and syn- 
chronous rotation of both components. We find a = 47 Rq, i = 60.3 
degree, stellar potentials 3.95 and 4.15 for the gainer and donor, re- 
spectively. This solution implies Mi = 27.7 Mq, M2 = 20.8 Mq, 
Ri = 15 Rq, i?2 = 1 1.8 Rq. These photometrically derived masses 
and orbital separation are consistent with the spectroscopic values 
given in Table 4 for an inclination angle of ~ 62 degree. 

The equivalent widths of lines of both stellar components are 
highly variable and roughly anticorrelated. The sum of both does 
vary as well, but with much lower amplitude. Unfortunately, with 
the available data it cannot be decided whether Pi or P2 does sort 
the data better. Similarly, the line widths are variable (see Fig. |4]l, 
with no clear preference for one of the two periods. 



5 DISCUSSION 



http://www.asu.cas.cz/ had/korel.html 



|Ofir| ( |2008[ l originally suggested a quadruple star in 
(Aa-l-Ab)-l-(Ba-l-Bb) configuration to explain the light curve of 
BI 108. Similarly, the light curve of CzeV343 was modeled as a 
double eclipsing binary with periods near a 3:2 ratio jCagas &] 
|Pejcha|[2012[ >. These authors speculate that this resonance might 
be the result of Kozai cycles and tidal friction during the evolution 
of a quadruple system, but recognize the lack of theoretical work 
to confirm this conjecture and that a simple coincidence cannot be 
discarded. In the case of BI 108, the proximity to perfect resonance 
is 0.008% (compared with 0.1% for CzeV343) but the circularity 
of the orbit is in principle inconsistent with Kozai cycles; they tend 
to increase the ellipticity of central binaries in multiple systems. 

For BI 108 the absence of color variations would require the 
Ba-l-Bb system has equal components. To make Ba-l-Bb as dim as 
possible (in order to hide eventual spectral features), one has to 
demand that the eclipse in the light-curve is total, such that the 
entire Ba-nBb subsystem contributes about 20% to the total flux 
as indicated by the photometric variability. Under such assump- 
tions the individual magnitudes would be 14.3 for Aa and Ab each, 
and 15.8 for Ba and Bb each, i.e. the subsystem A providing about 
80% of the total light in the /-band light curve. Rapid rotation (be- 
yond corotation of the order of 200 km s^^) for both components 
of Ba-l-Bb would still be required in order to hide spectral features 
from this pair from detection in the high S/N data. Taking all the 
above together, such a quadruple system cannot be ruled out and in 
spite of the lack of theoretical understanding, this will still turn out 
to be the hypothesis matching the requirements of Occam's razor. 

Another possible multiple system would consist of three stars 
in non-hierarchical orbits A-l-B-l-C, where Pi is the result of eclipses 
of the brighter component A by both B and C. Indeed some stable 
solutions exist for the unrestricted three body problem, for instance 
those found by Lagrange where three bodies orbit at the vertices of 
a rotating equilateral triangle. However, this would not only require 
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Figure 4. The spectral variability of Sim 4553, Hel4026 and H/S, sorted with the orbital period, from bottom to top. Offset is proportional to phase and 
according to values given in Table 1. The RVs as measured in H/9 are shown as filled disks for star#l and diamonds for star#2, respectively. 



a parameter fine-tuning beyond belief, but probably would still not 
be stable in a real stellar system undergoing stellar evolution. 



If there would be only two stars, they would show variabil- 
ity somehow resonantly synchronized with the orbit. Whatever this 
mechanism is, to achieve it, the two components must exert force 
on each other. This could be done through magnetic fields modify- 
ing the chemical element distribution on the stellar surfaces, hence 
producing line and continuum flux variability ( Krticka et al. 2007 1. 
Historically, it has been assumed that magnetic fields in OB stars 
are very rare if not absent. However, according to |Petit|f201 1^ , "... 
the increasing discoveries of fields in early B-type stars on the main 
sequence and pre-main sequence and in both young and evolved 
O-type stars show convincingly that fossil fields can and do ex- 
ist in stars with masses as large as 45 solar masses". However, in 
BI 108 magnetic fields can be excluded with high probability, due 
to the missing signatures of chemical inhomogeneities in the spec- 
tral lines. 



It has been shown that excitation of resonant non-radial g- 
mode pulsations should be a common phenomenon during tidal 
evolution of eccentric double main sequence binaries and it is also 
an effective mechanism for orbital decay ( |Witte & Savonije|1999| 
|2001 ). Although this is a promising mechanism for explaining the 
Pi variability in principle, it is difficult to imagine how it could 
generate the amplitude observed in the Pi light curve. This would 
require strongly perturbed stars, the visible stellar surface in the 
system varying by as much as 20 % between maximum and min- 
imum. Thus, other than for the double binary hypothesis, there is 
a theoretical framework for the variability, however the observed 
strength of the variation goes far beyond what has been explored in 
this framework. For this reason we consider it a less likely option 
than the double binary hypothesis, but note that neither can it be 
firmly ruled out. 



6 CONCLUSIONS 

We have presented the first spectroscopic study of the LMC system 
BI 108 showing an unique 3:2 resonant photometric variability. The 
immediate goal of the study was to determine which of the two 
periods (if not both) was due to an eclipsing binary. From the shape 
of the light-curves, it was originally expected that Pi would turn 
out to be the binary period. Surprisingly the radial velocities point 
to a SB2E system orbiting with orbital period P2 instead. 

The mechanism responsible for Pi remains unknown, but with 
the present data we exclude a system with three components and 
magnetic interaction between two binary components. The 3:2 res- 
onant tidal bulge would require amplitudes beyond what can realis- 
tically be expected from such a mechanism. The remaining sugges- 
tion, a double binary system seen by chance (or eventually by still 
unknown dynamical reasons) with a period ratio near 3:2, remains 
as the most plausible hypothesis. 
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